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ABSTRACT 

Using the extended J, H and K magnitudes provided by the 2MASS data archive, we consider the 
position of brightest cluster galaxies (BCGs) in the observed relations between inferred supermassive 
black hole (SMBH) mass and the host galaxy properties, as well as their position in the stellar velocity 
dispersion and luminosity (er* — L) relation, compared to E and SO galaxies. We find that SMBH 
masses (M B ) derived from near- infrared (NIR) magnitudes do not exceed ~3x 10 9 M Q and that these 
masses agree well with the predictions made from cr*. In the NIR, there is no evidence that BCGs 
leave the cr* — L relation defined by less luminous galaxies. The higher SMBH masses predicted from 
V-band luminosities (A/. < 1O 1O 5 M0) are attributed to the presence of extended envelopes around 
the BCGs, however, this will need to be confirmed using deeper multiwavelength imaging. 
Subject headings: galaxies: elliptical — galaxies: evolution — galaxies: fundamental parameters — 
galaxies: photometry 



1. INTRODUCTION 

Whether the first galaxies were formed from initial 
large-scale condensations, or grew from an assembly of 
smaller bodies, still remains one of the most fundamen- 
tal questions in modern astrophysics. Studies of the 
most massive galaxies will provide important constraints 
on this. Similar considerations apply to supermassive 
black holes (SMBHs) as the masses of SMBHs corre- 
late w ith properties of the host bulge (|Ferrarese fc Fordl 
2005), i.e., the S MBH mass vs. bulge lumi nosity 
(M. - L) relation (jKormendv fc Richstonel 11995( 1. the 
SMBH mass vs. stellar velocity di s persion (M, — 
cr*) r elation ((Ferrarese fc Merritt] 120001 : iGebhardt et all 
2000) , and the SM BH mass vs. Sersic index relation 
(|Graham fc Driven l200l . 

As highly luminous massive galaxies found toward 
the centers of galaxy clusters, brightest cluster galaxies 
(BCGs) have received considerable interest. The sur- 
face brig htness profile s (SBPs) of BCGs are well fit by 
the same lSersicl (119631) law tha t describes less-luminous 
spheroids ((Graham et al.lfl99 6L apart from the outer- 
most regio ns which sometimes exhibit faint, extended 
envelopes ((Oemlerl 119761 : iBernardi et alj 120061 . hereafter 
B06). BCGs also appear to obey the same relations be- 
tween fitti ng parameters that c haracterize E/SO ga laxies 
generally (|Graham et al.|[l996h . lLauer et alj (|2006l here- 
after L06) noted that the M, — L relation, in the V-band 
(M. — Ly), predicts higher SMBH masses in BCGs than 
are predicted by the M, — cr* relation. B06 obtain similar 
results; the slope in the size-luminosity relation is found 
to be steeper in BCGs when compared to the bulk of 
E/SO's, and the cr* vs. luminosity (cr* — Lr) relation is 
seen to flatten for the brightest galaxies. 

While the low-scatter M, — cr* relation is the pre- 
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ferred "secondary" SMBH mass estimation technique, 
compared to th e larg er-scatter M. — Ly relation, 
iMarconi fc Huntl (|2003l hereafter MH03) have shown 
that the scatter in the relations become similar if param- 
eters are derived in the near- infrared (NIR) . With this in 
mind we have conducted a study of BCGs based on the 
2MASS extended source catalog. We use the 219 L06 
galaxies of which ~ 30% are BCGs and the remainder 
are E/SOs. The L06 data include absolute V-band mag- 
nitudes (My) and, except in 51 value for cr, . We 
adopt the errors of 10% in My and cr*, as quoted by L06. 
We supplement the My data with the NIR data con- 
tained within the 2MASS extended source catalog. All 
magni tudes are corre c ted fo r galactic extinction accord- 
ing to Schlcge ~et al.l ([1998). Distances are all adjusted 
to a common scale, with Hq = 70 km s -1 Mpc - 1 , an d 
primarily taken from the survey of iTonrv et alj d200D. 
Remaining distances are taken from lLaine et al.l (2003) 
or from the Virgo in-f all corrected recessi onal velocities 
listed by Hyperleda 4 ((Paturel et al.ll2003H . In § [5] we 
evaluate the 2MASS photometry. In § [3] we present the 
results, which are discussed in §2J §[5] sums up. 

2. 2MASS PHOTOMETRY 

In this study we have used the 2MASS "total" magni- 
tudes (e.g., fc_m ext ) derived from SBP fitting extrapola- 
tion 5 , rather than aperture photometry which i nevitably 
under- estimates the total galaxy magnitudes ((Andreonl 
2002). Briefly, kjm cxt is estimated by numerical inte- 
gration of the Sersic law, fitted between r > 7-10" (to 
avoid the point spread function) and the maximum ra- 
dius (r max ) of the SBP with a signal to noise greater 
than two. Assuming circular isophotes and r max = 20- 
80" (see Figure [lj , this corresponds to a SBP limit from 
2.9 to 3.7 magnitudes below the 2MASS 3cr limit (20.09, 
19.34 and 18.55 in J, H and K). The best-fitting Sersic 
law is then integrated up to r ext (the galaxy "total" ra- 
dius). Typically, r cxt is ~ 2-5 times the radius of the 20 



4 http://leda.univ-lyonl.fr 

5 For more information "see sections 2.3a and 4.5e of 
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Fig. 1. — Simulations of 2MASS photometry. The distribution 
of simulated galaxy properties, total K-band magnitude and r ex t, 
(black circles) with respect to observed BCGs (red circles). 



80^ 

^60H 

^50^ 
us : 
OS 40 H 

1 30^ 
20 ^ 
W 10 1- 



«b 



10 20 30 40 



50 

2 , 



60 70 



SO 



R-band 23 mag/arcsec Radii (kpc) 



Fig. 3. — Comparing the radii of V-band and NIR, magnitudes. 
The solid line is one-to-one. Open and closed circles show the 
relation between r r 23 and r^o an d fext respectively. 
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Fig. 2. — Simulations of 2MASS photometry. The difference 
between "true" and estimated magnitudes from ~ 3500 simulations 
is shown. The dashed line marks the median of the distribution. 



mag/arcsec 2 isophote (e.g., rk2o) where 5rk20 is imposed 
as a strict upper limit. 

To quantify the importance of undetected light at large 
radii, we have carried out extensive simulations of BCGs 
and performed 2MASS analogous photometry. Over 
3500 BCGs with My and z randomly distributed in the 
observed ranges (-24.5 - -22 and 0.015 - 0.050, respec- 
tively) were generated. For each BCG the effective radius 
was derived from the r c — My relation provided by L06, 
and the Sersic index, n, was deriv ed from the n — r e re - 
lation estimated from Figure 11 of iGraham et~aT1 (|l996h 
(0.2 dex scatter was included for both). In all cases we 
imposed 1 < n < 15. The observed total K magnitude, 
Mr, was then derived using V-K=3.6. The model images 
of each BCG (characteriz ed by wki ? * e and n) were then 
generated using GALFIT ()Peng et all l2002f) to take into 
account the 2MASS pixel sizes (1") and spatial resolu- 
tion (FWHM - 3"). We added 2MASS typical noise (K 
= 19.74 mag/arcsec 2 rms) and derived SBPs that were 
fitted with a Sersic law to estimate r cxt and fc_m oxt . To 
reproduce the magnitude-size relation (see Figure [1} we 
used an average r cxt of 2.4rfc 2 o- Figures Q] and [5] show how 
the simulated magnitudes compare to the actual BCG 
magnitudes and the "true" simulated input magnitudes. 
The average offset is -0.5 mags (50 th percentile); 2MASS 
mildly underestimates Mk in BCGs. Simulations using 
bluer colors (V-K=3.3) give very similar result. 

In § |3] we will directly compare the results from ex- 
trapolating photometry estimated in the V-band and 
NIR. As the V-band photometry may have been gained 
from different radii (thereby gathering a different amount 
of light) we now briefly compare the relative sizes of 
the regions from which the photometry was derived. 
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Fig. 4. — The relationship between <r„ and luminosity in the 
V-band (a) and K-band (b). Open black circles mark E and SOs, 
filled red circles show BCGs. The solid li nes show the fit to ju st 
the E/SOs, the dotted lines mark the fit of lBernardi et al.l (2006). 



The L06 V - band BCG luminosities were provided by 
lLaine et"aH ((2003). who in turn es timated the total 
magni tudes from SBPs presen ted by iPostman fc LaueJ 
(|1995h . IGraham et all |l996f ) performed Sersic fits to 
these SBPs using a limiting R-band surface brightness. 
Therefore, assuming a color of R-K=3.0, we can directly 
compare the radii of r oxt and r^o (the radius of the K=20 
mag/arcsec 2 isophote) with the R-band surface bright- 
ness at 23.0 mag/a r csec 2 (r r 2a). By calculating r r2 3 from 
the IGraham et ail (|1996D BCG Sersic fits, we find that 
fk20 agrees well with the r r 23 radii (Figure [3]). However, 
r ext is always larger than r r 23- Therefore, the 2MASS 
magnitudes used in this study include light from at least 
as extended a region as the V-band magnitudes used by 
L06. 

As an additional check, we can compare published 
values of Mk to those derived by MH03 from two- 
dimensional fitting to 2MASS profiles out to infinity. We 
find a systematic offset of ~ 0.4 mags, consistent with 
our simulations. 



NIR Luminosities in BCGs 
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Fig. 5. — V-band vs. NIR Mm estimates. The solid line marks 
a one-to-one relation. Open circles are E/SO's, closed circles are 
BCGs. Blue, green and red colors refer to the J, H and K bands 
respectively. The dotted line marks the fit to all the K-band data. 

3. BCGS IN THE NIR 

Figure [4] compares the cr* — Ly and cr* — Lk relations. 
In both cases we plot the best-fit relation defined by the 
E/SO galaxies (solid line) as well as the cr* — Lr fit given 
by B06 (their Figure 6 with colors of V-R=0.6 and R- 
K=3.0). We find shallower slopes for the E/SO popula- 
tion consistent with B06. Figure BJa) demonstrates that 
the "bending" of the cr* — Ly relation, as noted by B06, 
is also seen in the L06 sample; BCGs fall above the re- 
lation defined by the E/SOs. In the NIR (Figured)) the 
BCGs do not appear to define a separate population; in- 
stead their distribution is indistinguishable from that of 
the E/SOs. The average offset of BCGs from the E/SO 
relation is 1.20 mags in My and 0.48 mags in Mk. 

Figure \5\ presents the relationship between M, esti- 
mated from the V-band (hereafter M.(V)), using the re- 
lation as defined by L06 (their Equation 4), and M. esti- 
mated from the NIR data (hereafter M.( J, H, K)) using 
the relations defined by the MH03 sample. The upper 
limit for Af.(V) is 10 la5 M Q . Below M. « 10 8 5 M Q the 
agreement between all bands is good. Above 1O 8 ' 5 M , 
the NIR data predict significantly lower SMBH masses, 
with none exceeding 10 9 4 M Q . The fit to the K-band 
relation is shown as a dotted line and has a slope of 
0.62 ± 0.02 (0.74 ± 0.02 for E/SOs). For estimates of 
M. (J, H, K) we do not use the exact fits presented by 
MH03 (e.g. their Table 2) as they were derived using 
a bisector method (e.g., assuming scatter in both M, 
and £j,h,k)' Instead, M,(J, H, K) has been determined 
from a single Y|X fit to M, and L (taking into account 
errors from both axis) because in estimating M. from L 
one assumes that all the scatter is in only one variable. 
The form of this relation, in the K-band, is given by 
logM. =8.22±0.07+(1.06±0.11)(logL K -10.9). The 
results, in this case, are insignificant from the original 
MH03 fits (cr = 10 04 M Q ). 

Figure[6]shows how photometric M, est imates compare 
to tho se from the M. — cr* relation. The iTremaine et al.l 
( 2002) expression is used to derive M. from cr* , hereafter 
M.fcr*). In Figure [6ja) the M, — Ly relation (above 
10 8 °M Q ) predicts SMBH masses greater than those ex- 
pected from the M m — cr, relation. However, in the NIR, 
this observation is not made; both predictions are consis- 
tent. The M. - L v relation implies M, < 2.5 x 1O 1O M , 
whereas the NIR produces M. < 2.8 x 10 9 M Q . The scat- 
ter in the M.(J, H, K)—M,(a*) relations are significantly 
less than the M,(V) — Al.(cr») relation. 



The effect of the results from the § [2] simulations can 
be seen by artificially and randomly introducing the dis- 
tribution of A mag to the K-band magnitudes. A com- 
parison of M, (cr*) and the adjusted M,(K) is shown in 
Figure [7] The slope of the best fit relation (dotted line) 
is 0.85 ± 0.05 and no values of M.(K) exceed 10 9 7 M Q . 

4. DISCUSSION 

It is evident from both Figures[3]and[H]that the disper- 
sion of BCGs in the NIR is considerably less than in the 
V-band. This is also the case for the E/SOs. Even if we 
are underestimating the NIR luminosities, it would re- 
quire a very fortunate coincidence to have such a strong 
agreement between M, (cr») and M,(J, H, K) across the 
entire mass function. It then follows that BCGs are not 
"special" when viewed at NIR wavelengths. BCGs follow 
the same cr* — Lk distribution as defined by less luminous 
spheroids, and comparable masses are predicted for the 
SMBHs in BCGs based either on velocity dispersions or 
on total magnitudes. 

We have demonstrated, through extensive simulations, 
that the 2MASS magnitudes used in this study are ro- 
bust to within 0.5 mags for BCGs. This is consistent 
with the offset between the 2MASS photometry and that 
of MH03, who take the total magnitude components of 
Sersic fits to two-dimensional photometric models. This 
offset is expected as the 2MASS total magnitude integra- 
tions stop at r ext . For the M, — Lj^k relations to predict 
a population of 10 10 Af© SMBHs, i.e., for the BCGs to 
fall on the one-to-one relation in Figure El 2MASS pho- 
tometry would have to be increased, on average, by 1.65 
mags. This offset is inconsistent with the 2MASS under- 
estimates. The offset would also need to vary with M, 
and have no effect at masses below 1O 8 ' 5 M where the 
results between M, — Ly and M, — -L^h.k are consis- 
tent. Conversely, the V-band photometry must be over- 
estimated by 2.17 mags, on average, for the upper limit 
of M, to be similar to that predicted from the NIR. 

Why then do BCGs show an excess of V-band light 
over their less massive cousins? The issue could be re- 
solved by considering the faint extended luminou s halos 
know n to surround the most massive galaxies (|Oemlerl 
1976). The signature of such halos is an inflection in 
the SBPs at large r. In the NIR, the total magnitudes 
could be missing the contributions from faint halos or, 
alternatively, V-band luminosities could be interpreted 
as over-estimates due to the spurious inclusion of halo 
light. Since extended halos predominantly sit in the over- 
all cluster potential, they are unlikely to be related dy- 
namically to the central regions from which the BCG 
cr» is typically measured. Photometry from further to- 
ward the blue end of the spectrum may be deep enough 
to include a significant contribution from these extended 
halos, leading to an increase in the estimation of BCGs 
luminosities and a turn-over in M. estimates above a 
certain threshold. 

While we have shown that BCGs are not "special" in 
terms of their cr* — Lk distribution, nevertheless, their 
extreme luminosities and unique locations close to the 
centers of galaxy clusters suggest special formation pro- 
cesses. It has long been argued that the extended en- 
velopes of BCGs are debris from tidally-stripped galax- 
ies, and hence that they are associated more closely with 
the overall cluster potential well than with any single 
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Fig. 6. — Comparing photometric M, estimates with M, — <t» estimates. V, J, H and K masses (a, b, c and d respectively) are compared 
to M. — rj„ estimates. Open black circles are E/SOs, BCGs are filled red circles. In all cases the solid line represents a one-to-one relation 
and the dotted line the best-fit relation. 
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Fig. 7. — Comparing redistributed M,(K) estimates with M. — 
o t estimates. Black circles are E/SOs, BCGs are red circles. The 
solid line represents a one-to-one relation and the dotted line the 
best fit relation. 

galaxy (|Merrittlll984a| L The envelopes m ay also consis t 
in part of stars formed in cooling flows (Fabian 1994). 
The presence of multiple nuclei in some BCGs a rgues in 
favor o f these galaxies not being fully relaxed (|Merrittl 
1984b). A photometrically complete, high-resolution 
imaging survey of BCGs would be able to provide a 
framework for a more quantitative analysis of these fun- 
damentally important objects. 

5. CONCLUSIONS 

Brightest cluster galaxies offer the chance to study the 
pinnacle of galaxy evolution. They also give us the op- 
portunity to study the top of the SMBH food chain by 
using the observed relations between M. and the prop- 
erties of the surrounding host galaxy. We have shown 
that NIR luminosities, combined with previously estab- 
lished scaling relations (MH03), imply a maximum mass 
of ~ 3 x 10 9 M©. This is consistent with the most 



massive SMBH directl y modeled at the center of M87 
(iMacchetto et al.lll997fl and with the direct M, estimates 
of IDalla Bonta et all (12006) in 3 BCGs. We also find 
that, across all values, SMBH masses predicted using 
NIR magnitudes are consistent with masses predicted 
from cr*. In addition, we have shown that BCGs follow 
the same distribution, as defined by E/S0 galaxies, in the 
a„ —Lk relation. If confirmed by a deep multiwavelength 
study, these findings could have important implications 
for the nature of the SMBH mass function, and, as in 
the past, would show that NIR data are to be preferred 
when estimating M. (MH03). While BCGs are likely 
not special - in the sense of hosting hyper-massive black 
holes or by defining a distinct population in the 0% — Lk 
plane - they may be interesting by virtue of being sur- 
rounded by extended faint halos. The unique local BCG 
environment, deep within a cluster potential, could be 
the generator of these halos, which may be populated by 
younger stars tidally stripped from other cluster mem- 
bers or that are the results of intra-cluster gas accretion 
or other recent merger events. 
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